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1 .  Introduction 

Fiber  Bragg  Grating  (FBG)  sensors  are  an  increasingly  important  emerging  technology  in 
the  area  of  intelligent  structural  health  monitoring  (SHM)  of  civil,  mechanical,  naval  and 
aerospace  structures.1'9  A  large  number  of  FBG  sensors  can  be  easily  written  in  a  single 
fiber.  For  SHM  applications,  FBG  sensors  are  typically  used  to  monitor  static  or 
dynamic  strains.  The  strains  cause  spectral  shifts  in  the  FBG  sensor  transmittivity 
(equivalently  reflectivity).  The  spectral  shifts  can  be  monitored  in  several  ways  with  an 
appropriate  demodulator.  For  spectrally-encoded  FBG  sensors,  extant  demodulation 
schemes  can  be  classified  into  three  categories,  namely:  scanning  type,  spectrometry- 
based  and  interferometry-based.  Scanning  type  techniques  include  Fabry-Perot  scanning 
filters,10  acousto-optic  tunable  filters,11  and  tunable  laser  sources.12'13  All  these  scanning 
type  techniques  suffer  from  the  fact  that  at  any  instant  only  one  FBG  sensor  can  be 
interrogated.  Such  approaches  are  not  applicable  if  all  the  sensors  have  to  be  interrogated 
simultaneously  for  the  purpose  of  monitoring  impact  signals  and  acoustic  emissions.14 
Spectrometric  methods15  suffer  from  low  sensitivity  and  are  not  suitable  for  dynamic 
measurements  if  several  sensors  have  to  be  active  at  all  times.  Interferometric  methods 
such  as  the  Mach-Zehnder  interferometer16  are  ideally-suited  to  monitor  dynamic  strains; 
however  they  require  electronic  feedback  to  actively  compensate  for  any  quasistatic  drift 
to  maintain  the  interferometer  at  quadrature.  This  makes  the  cost  of  multiplexing  high 
since  each  sensor  requires  its  own  feedback  system.  A  cost-effective  and  parallel 
demodulation  scheme  for  arrays  of  FBG  sensors  is  therefore  necessary. 

In  work  over  the  past  two  years,  we  have  demonstrated  a  novel  two-wave  mixing 
(TWM)  wavelength  demodulation  scheme  for  FBG  sensors  that  has  the  ability  to 
compensate  for  quasistatic  drifts  without  the  need  for  active  stabilization.  This  method  is 
suitable  for  measuring  dynamic  and  transient  strains  such  as  those  induced  by  vibrations, 
impact,  ultrasound,  or  acoustic  emission.  Furthermore,  this  scheme  can  be  readily 
multiplexed  by  using  wavelength  division  multiplexing.  Thus  far,  a  four-channel  TWM 
device  has  been  demonstrated  to  demodulate  dynamic  strain  signals  from  four  FBG 
sensors  simultaneously. 

In  Section  II,  the  details  of  accomplishments  to  date  are  given.  In  section  III,  a 
plan  of  work  for  continued  research  is  outlined.  The  aim  of  the  continuation  phase  is  to 
scale  up  the  system  developed  so  far  and  to  demonstrate  application  on  a  real  structure 
with  a  view  to  transitioning  the  technology. 


n.  Work  to  Date: 

2.1  Adaptive  Two-wave  Mixing  Photorefractive  two-wave  mixing  is  essentially  a  dynamic 

holographic  process  where  two  coherent  optical  beams, 
typically  one  stronger  pump  beam  and  one  weaker  signal  beam,  interact  within  a 
photorefractive  crystal  (PRC).17'18  The  process  can  be  briefly  summarized  as:  (1)  creation 
of  intensity  gratings  due  to  coherent  stationary  interference  of  the  interacting  beams 
leading  to  (2)  non-uniform  photo-excitation  of  charges  in  the  PRC  which  then 
diffuse/drift  to  create  (3)  a  space-charge  field  within  the  PRC,  which  in  turn  creates  (4)  a 
refractive  index  grating  via  the  electro-optical  effect,  and  which  causes  (5)  diffraction  of 
the  interacting  beams.  This  process  of  dynamic  coupling  is  known  as  wave  mixing,  a 
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phenomenon  that  has  been  exploited  over  the  years  for  a  variety  of  applications17'18 
including  in  optical  interferometry.19'24 

Figure  1  shows  a  schematic  of  a  simple  TWM  interferometer  in  the  so-called 
direct  detection  configuration21  used  in  this  paper.  The  signal  beam  and  pump  beam 
interfere  with  each  other  inside  the  PRC  and  create  an  index  grating  in  the  PRC.  A  part  of 
the  pump  beam  is  diffracted  by  the  created  index  grating  into  the  direction  of  the 
transmitted  signal  beam.  The  diffracted  pump  beam  is  wavefront-matched  with  the 
quasistatic  signal  beam.  Subsequent  to  the  PRC,  the  diffracted  pump  beam  and  the 
transmitted  signal  beam  interfere  with  each  other  to  demodulate  any  phase  difference 
between  them.  The  TWM  interferometer  is  adaptive  because  the  created  index  grating  is 
a  dynamic  grating  and  the  crystal  can  adapt  to  any  phase  shift  that  is  slower  than  the  PRC 
response  time  by  forming  a  new  index  grating.  The  diffracted  pump  beam  will  therefore 
track  any  quasistatic  changes  in  the  signal  beam  phase,  resulting  in  no  net  quasistatic 
phase  difference  between  the  diffracted  pump  and  transmitted  signal  beams.  Only 
dynamic  phase  changes  faster  than  the  photorefractive  response  time  of  the  PRC  will  not 
be  present  in  the  diffracted  pump  beam.  This  will  result  in  a  net  phase  difference  between 
the  diffracted  pump  and  the  transmitted  signal  beams  subsequent  to  the  PRC,  which  can 
then  be  picked  up  by  interference  of  these  two  beams.  By  application  of  an  external  DC 
field  to  the  PRC  (the  so-called  drift  regime),  the  interference  pattern  can  be  kept  nearly  in 
phase  with  the  created  index  grating.  5  in  the  drift  regime,  by  adjustment  of  the  applied 
field  and  two  beam  mixing  angle,  the  TWM  interferometer  can  be  made  to  work  at  near 
quadrature  (that  is,  the  diffracted  pump  beam  is  nearly  jr/2  phase  shifted  with  respect  to 
the  transmitted  signal  beam)  for  sensitive  and  linear  interferometric  detection.  The 
adaptivity  of  the  TWM  process  in  the  PRC  implies  that  the  system  will  remain  at 
quadrature  without  any  need  for  active  external  stabilization  despite  quasistatic  drifts  that 
may  occur.  The  PRC  can  therefore  be  regarded  as  an  adaptive  high-pass  filter  that 
selectively  monitors  any  high  frequency  changes  in  the  mixing  beams.  This  is  ideal  for 
monitoring  small  dynamic  strain  signals  (such  as  those  induced  by  impact  or  acoustic 
emissions)  that  may  ride  on  top  of  large  quasistatic  drifts  (caused  by  static  stressing  or 
temperature).  It  should  also  be  noted  that  the  TWM  interferometer  is  readily 
multiplexable  using  any  combination  of  angular  multiplexing,  spatial  multiplexing  or 

This  enables  multiple  interferometers  to  be  configured  in  one 


wavelength  multiplexing, 
simple  device.26'2 


Figure  1  Schematic  of  the  TWM  direct  detection  setup 
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Following  the  analysis  of  Delaye  et  al21,  the  amplitude  of  the  transmitted  signal  beam 
subsequent  to  the  PRC  can  be  denoted  as: 

Es  =  E0  exp[z'(2>(/)]  (1) 

Here  <j>(t)  is  the  dynamic  phase-shift  which  is  the  signal  of  interest.  E0  is  the  complex 
amplitude  of  the  signal  beam  and  may  include  an  additional  phase  term  <p„  to  account  for 
quasistatic  phase  shifts  due  to  noise.  The  amplitude  of  the  diffracted  pump  beam  EdP  can 
then  be  written  as  a  function  of  the  TWM  gain  as  follows:21 

Edp  =£°{exp[^]-1}’  (2) 

in  which  y  =  y  +  iy"  is  the  TWM  complex  gain  and  1  is  the  crystal  length  in  the  beam 
propagation  direction.  It  is  the  interference  between  the  diffracted  pump  edP  and  the 
transmitted  signal  beam  es  that  demodulates  the  spectral  shift  .  In  the  limit  of  small 
dynamic  phase  shifts  (this  is  usually  valid  for  the  dynamic  strains  induced  by  acoustic 
emission  or  impact,  for  instance),  the  interference  signal  at  the  photodetector  can  be 
written  as: 

/  =  \E„\ 2 {e2^  +  2sin(y" L)er L  <p(t)}  (3) 

The  first  term  in  the  curly  brackets  represents  a  dc  level  and  contributes  to  photodetector 
shot  noise.  The  second  term  represents  the  signal  of  interest,  and  indicates  that  the 
interference  signal  varies  linearly  with  the  dynamic  phase  shift  (for  small  phase  shifts) 
between  the  signal  and  pump  beams.  It  is  clear  from  eq.  (3)  that  the  maximum  signal 
occurs  when  y”L  approaches  nil,  that  is  at  quadrature.  This  condition  can  be  readily 
approached  in  the  drift  regime  by  application  of  an  electric  field  to  the  PRC,  and  by 
appropriately  adjusting  the  TWM  parameters. 

2.2  Adaptive  Two-wave  Mixing  Wavelength  Demodulator  now  show  that  the  TWM 

interferometer  can  be 
configured  as  a  wavelength  demodulator  for  FBG  sensors.  The  light  reflected  from  an 
FBG  sensor  is  first  split  into  two  beams,  one  pump  and  one  signal  beam,  which  are  made 
to  travel  unequal  paths  prior  to  mixing  in  the  PRC.  Any  wavelength  shift  of  the  light 
reflected  from  the  FBG  sensor  then  results  in  an  equivalent  phase  shift  between  the  pump 
and  signal  beams  because  of  travel  over  unequal  optical  paths.  Therefore,  in  the  sense  of 
converting  a  wavelength  shift  of  the  FBG  sensor  signal  into  an  equivalent  phase  shift,  the 
underlying  principle  used  here  is  the  same  as  that  used  in  Mach  Zehnder  Interferometry 
and  other  unbalanced-path  interferometric  methods.16 

The  spectral  shift  induced  by  strain  and  temperature  in  a  Bragg-grating  sensor  can 
be  written  as:29 

2 

A/ln  neff 

— —  =  [l--^-tPi2-L'(Pn  +  Pi2)^z  +(«A  +aN) AT  ’  (4) 

AB  2 

in  which  XB  is  the  center  wavelength  of  the  FBG  sensor,  AXB  is  the  wavelength  shift 
caused  by  strain  or  temperature,  neff  is  the  effective  refractive  index  of  the  fiber,  ptj  are  the 
components  of  the  strain-optic  tensor,  v  is  Poisson’s  ratio,  sz  is  the  strain  along  the  fiber, 
cxa  is  the  thermal  expansion  coefficient  and  an  is  the  thermo-optic  coefficient  of  the  fiber. 
For  typical  Bragg-sensors  at  1550nm,  it  has  been  estimated  that  one  micro-strain  (lps 
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=10'6)  will  lead  to  a  ~1.2  pm  change  in  wavelength,  and  a  1°C  change  in  temperature  will 
lead  to  about  13pm  change  in  wavelength.29 

As  mentioned  earlier,  a  spectral  shift  AABof  the  signal  and  pump  beams  is 
effectively  converted  into  a  relative  phase  shift  between  the  beams  (due  to  travel  through 
unbalanced  optical  paths),  and  this  is  given  by:29 

(5) 

A B 

where  d  is  the  optical  path  difference  (OPD),  As  is  the  center  wavelength  of  the  FBG 
sensor,  Mb  is  the  wavelength  shift  of  the  light  reflected  by  the  FBG.  From  Eq(5),  it 
appears  that  the  greater  the  OPD,  the  larger  the  equivalent  phase  shift,  and  therefore  the 
stronger  the  interference  signal  should  be.  However,  typically  broadband  light  sources 
are  used  to  illuminate  the  FBG  sensors,  and  the  FBG  reflection  spectrum  typically  has  a 
finite  line  width  on  the  order  of  0.1-0.4nm.  This  implies  that  coherence  of  the  two 
interfering  beams  needs  to  be  taken  into  account  both  in  the  photorefractive  grating 
creating  process  and  in  the  subsequent  interference  between  the  transmitted  signal  and  the 
diffracted  pump  beams.  The  fringe  visibility  due  to  interference  of  two  beams  of  finite 
spectral  width  A k  is  given  by: 


V  =  x  exp(. 
r  + 1 


mV. 

16  In  2  ' 


where  r  is  the  intensity  ratio  of  the  two  beams,  and  d  is  the  OPD.  Incorporating  the 
degradation  in  fringe  visibility  due  to  low  coherence  in  the  TWM  analysis  leads  to  a 
modified  expression  for  Eqn(3).  The  TWM  spectral-shift  demodulated  interference 
signal  (neglecting  the  dc-component  for  simplicity)  can  be  shown  to  be  proportional  to: 


S  oc  exp(/Z)sin(y"i)  exp(- 


Ak2d2  d AAb 
'  161n2  4 


The  effect  of  decreased  fringe  visibility  on  the  photorefractive  grating  formation  process 
is  not  a  significant  factor,  since  the  preferred  mode  of  two- wave  mixing  in  PRCs  is  in  the 
low  intensity  modulation  regime  (this  avoids  the  creation  of  higher  order  index  gratings 
leading  to  multiple  diffracted  beams).  Furthermore,  the  system  is  operated  at  near 
quadrature,  y'L&nl  2,  and  since  the  TWM  energy  gain  (the  real  part  of  the  complex 
gain)  i$  typically  very  small,  the  variations  in  the  first  two  terms  with  the  OPD  are  not 
significant  (as  will  be  demonstrated  in  the  experimental  section).  The  dependence  of  the 
interference  signal  on  the  OPD  is  therefore  predominantly  given  by: 

~  0  f  Ak2d2\dAAB  - 

S  oc  exp( - >  — r-2-  (8) 

I  161n2  XL 


The  wavelength  demodulation  signal  is  therefore  a  strong  function  of  the  OPD  d  with  the 
exponential  decay  arising  from  loss  of  coherence,  and  the  linear  increase  arising  from  the 
increase  with  OPD  in  the  phase-shift  due  to  spectral  change.  Figure  2  is  a  plot  of  the 
signal  amplitude  versus  the  OPD  for  different  line-width  FBG  sensors.  (In  the  calculation 
of  Figure  2,  L  was  set  equal  to  1cm  and  a  value  of  0.3  cm'1  for  was  used  for  the  real  part 
of  the  TWM  complex  gain).  As  indicated  in  Figure  2,  when  the  OPD  equals  zero,  there  is 
no  wavelength-demodulated  signal  detected,  and  as  the  OPD  increases  the  signal 
amplitude  increases  to  a  maximum  beyond  which  the  signal  starts  to  drop  due  to 
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decreasing  fringe  visibility.  For  each  given  line  width  of  the  FBG  sensor,  there  exists  an 
optimum  value  of  the  OPD  that  maximizes  the  wavelength-demodulated  signal.  It  is  also 
clear  from  Figure  2  that  the  narrower  the  line  width,  the  larger  the  optimum  OPD  and 
therefore  the  larger  the  demodulated  signal.  However,  it  is  not  always  better  to  use  a 
narrower  line-width  FBG  sensor.  First,  a  narrower  line-width  FBG  sensor  is  usually 
longer  in  length,29  which  decreases  the  highest  frequency  to  which  the  FBG  can 
respond.30  Secondly,  a  larger  OPD  decreases  the  dynamic  range  that  the  FBG  sensor  can 
measure. 


Figure  2  Wavelength  demodulated  signal  amplitude  as  a  function  of  the  optical  path 
difference  (OPD)  for  different  line-width  FBG  sensors. 


2.3  Experimental  Configuration  The  experimental  configuration  is  shown  in  Figure  3. 

The  FBG  sensor  is  illuminated  by  a  broadband  amplified 
spontaneous  emission  (ASE)  source  in  the  C-band  (1530nm  to  1570nm),  and  the 
reflected  light  is  coupled  by  a  circulator  into  an  Erbium  doped  fiber  amplifier  (EDFA) 
that  works  in  saturation  mode  (output  500m W).  The  amplified  light  is  split  using  a  1x2 
coupler  (splitting  ratio  95/5)  into  pump  (95%)  and  signal  (5%)  beams  that  travel 
unbalanced  optical  paths  to  the  InP:Fe  photorefractive  crystal.  The  InP:Fe  crystal  is 
oriented  for  two-wave  mixing  in  the  direct  detection  configuration,21  in  which  both 
beams  enter  the  ciystal  by  the  [1 10]  face.  An  external  DC  field  is  applied  along  <00 1> 
direction.  In  order  to  apply  a  continous  DC  field  across  the  InP:Fe  ciystal,  a  peltier  cooler 
is  used  to  prevent  electrical  breakdown  due  to  crystal  overheating.  (In  some  of  the  earlier 
experiments  reported  here,  the  peltier  cooler  had  not  been  integrated  into  the  set  up,  and 
in  these  cases  the  DC-field  was  applied  intermittently  to  avoid  possible  damage  to  the 
crystal  due  to  heating.)  Two  half-wave  plates  (HWP)  are  used  prior  to  the  PRC  to  rotate 
the  pump  and  signal  beam  polarization  to  be  s-polarized  (along  the  <1 10>  direction).  Due 
to  interference  between  the  signal  and  pump  beams,  a  refractive  index  grating  forms  in 
the  PRC  along  the  <00 1>  direction.  Under  the  applied  external  DC  field  (drift  regime), 
the  index  grating  is  nearly  in  phase  with  respect  to  the  interference  pattern,  and  the 
diffracted  pump  and  the  transmitted  signal  beams  are  kept  at  near  quadrature.  This  is 
optimal  for  demodulation  of  vety  small  phase/wavelength  changes. 
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Figure  3  Experimental  configuration  of  the  FBG  sensor  and  the  two-wave  mixing  (TWM) 
wavelength  demodulator. 

2.4  Wavelength  Demodulation  In  a  first  set  of  experiments,  an  FBG  sensor  centered  at 

1552nm  with  a  line-width  of  0.1  nm,  length  of  10mm,  and 
reflectivity  of  50%  was  glued  onto  a  PZT  stretcher  that  was  used  to  induce  known 
dynamic  strains.  The  light  reflected  from  the  FBG  sensor  undergoes  spectral  shift  due  to 
strain-induced  changes  in  the  Bragg-reflectivity.  The  sensor  by  itself  is  sensitive  to  both 
quasistatic  and  dynamic  strains,  and  is  also  subject  to  temperature  drift.  Flowever,  as 
explained  in  the  previous  section,  the  TWM  demodulator  system  will  automatically 
compensate  for  quasistatic  drifts  and  track  only  the  dynamic  strains. 


0  2  4  0 

Time(ms) 

Figure  4  Wavelength  demodulated  signal  at  different  values  of  OPD.  An  intermittent  DC 
field  is  applied  to  the  PRC  starting  from  1ms  to  6ms  and  a  10  kHz  10  pe  strain  is  applied 
to  the  FBG  from  2ms  to  6ms. 
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In  order  to  demonstrate  wavelength  demodulation,  we  applied  a  10  kHz,  10  ps 
strain  onto  the  FBG  sensor  and  measured  the  wavelength  demodulated  signal  amplitude 
at  different  values  of  the  optical  path  difference.  As  shown  in  Figure  4,  an  intermittent 
DC  field  is  applied  from  1ms  to  6ms  with  respect  to  a  reference  trigger,  and  the 
photorefractive  grating  initially  builds  up.  The  dynamic  strain  is  applied  as  a  toneburst 
starting  from  2ms  to  6ms.When  the  OPD  equals  to  zero,  although  the  TWM  energy  gain 
is  at  its  maximum,  there  is  no  detected  wavelength  demodulated  signal  because  there  is 
no  OPD  to  convert  the  wavelength  shift  into  phase  shift.  As  the  OPD  increases,  the 
wavelength  demodulated  signal  starts  to  appear. 


OPD(mm) 

Figure  5  Plot  of  wavelength  demodulated  signal  amplitude  and  TWM  gain  versus  OPD. 

Figure  5  is  a  plot  of  the  wavelength  demodulated  signal  amplitude  and  TWM 
energy  gain  versus  the  OPD.  The  optimum  OPD  that  maximizes  the  wavelength 
demodulated  signal  for  the  0.1  nm  line-width  FBG  sensor  is  found  to  be  8mm.  For  an 
OPD  of  8mm,  Eq.  (5)  indicates  a  wavelength-to-phase  shift  conversion  sensitivity  of 
about  *21  radians  per  nanometer  wavelength  shift  at  1550nm  wavelength.  This  translates 
to  0.0252  radian/pe.  Such  phase  shifts  are  readily  detectable  by  the  TWM  interferometer. 
Also  note  that  the  trend  of  the  signal  amplitude  curve  is  similar  to  that  of  the  theoretical 
curve  shown  in  Figure  2.  Also,  as  we  mentioned  in  the  theoretical  section,  the  TWM 
energy  gain  /is  indeed  a  small  number  and  experimentally  is  found  to  vary  from  0.47 
cm'1  to  0.1  cm'1  as  the  OPD  changes  from  0  to  12mm.  This  causes  the  first  exponential 
term  in  Eq.(7)  to  vary  from  1.6  to  1.1,  which  is  much  smaller  than  the  change  in  the 
second  exponential  term  which  decreases  from  1  to  0.03  as  the  OPD  increases  from  zero 
to  12mm. 

The  TWM  demodulator  output  when  operated  at  near  quadrature  is  expected  to  be 
a  linear  function  of  the  FBG  spectral-shift  (and  hence  strain).  This  is  shown  in  Figure  6 
where  the  demodulator  output  is  plotted  as  a  function  of  the  dynamic  (20  kHz)  strain 
amplitude.  The  minimum  detectable  strain  with  our  current  setup  was  measured  to  be 
~0.25pe,  corresponding  to  a  spectral  shift  of  0.3pm.  The  minimum  detectable  spectral 
shift,  which  is  limited  by  the  ASE  source  and  EDFA  intensity  noise,  can  be  further 
improved  by  using  balanced  photodetection  to  cancel  the  intensity  noise. 
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Figure  6  TWM  demodulator  output  as  a  function  of  dynamic  strain  amplitude. 

2.5  Adaptivity  to  Quasistatic  Drift  0ne  of  the  advantages  of  using  the  two-wave  mixing 

interferometer  as  a  wavelength  demodulator  is  its 
adaptivity  to  low  frequency  drift.  As  mentioned  earlier,  the  two-wave  mixing 
interferometer  is  automatically  adaptive  to  low  frequency  strain  or  temperature  drift  of 
the  FBG  sensor.  In  order  to  demonstrate  this  feature,  we  subjected  the  FBG  sensor  to 
quasistatic  strain  and  temperature  drift,  and  tested  the  system  response. 
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Figure  7  TWM  wavelength  demodulator  response  to  a  frequency  sweep  signal  from  10 
Hz  to  1.2  kHz. 
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Figure  8  (a)  Fourier  spectrum  of  the  applied  frequency  sweep  signal  from  10  Hz  to  1.2 
kHz.  (b)  Fourier  spectrum  of  the  response  of  the  TWM  wavelength  demodulator,  (c) 
Transfer  function  of  the  TWM  wavelength  demodulator.  The  cut-off  frequency  is  seen  to 
be  around  600  Hz  for  this  configuration. 

In  order  to  show  adaptivity  to  quasistatic  strain,  we  applied  a  frequency  sweep 
signal  from  10Hz  to  1.2  kHz  with  a  constant  magnitude  of  lOpe.  Figure  7  shows  the 
response  of  the  wavelength  demodulator  to  this  frequency  sweep  signal.  The  demodulator 
ignores  the  low  frequency  strain  applied  in  the  beginning  and  starts  to  respond  to 
frequencies  above  600Hz.  Figure  8  demonstrates  the  adaptivity  to  quasistatic  strain  more 
clearly  in  the  frequency  domain.  Figure  8  (a)  is  the  Fourier  spectrum  of  the  applied  sweep 
signal.  Figure  8  (b)  is  the  spectrum  of  the  system  response  to  the  sweep  signal,  and  Figure 
8  (c)  is  the  system  transfer  function  (modulus)  calculated  by  dividing  the  system  response 
spectrum  by  the  sweep  signal  spectrum. 

From  Figure  8,  it  is  clear  that  the  TWM  wavelength  demodulator  is  adaptive  to 
low  frequency  strains  and  acts  like  a  high  pass  filter  with  a  cut-off  frequency  of  600  Hz. 
The  cut-off  frequency  of  the  system  is  directly  related  to  the  response  time  of  the  InP:Fe 
PRC.  The  faster  the  PRC  (shorter  response  time),  the  higher  the  cut-off  frequency  will  be. 
Furthermore,  the  response  time  of  the  PRC  can  be  tuned  by  adjusting  the  input  optical 
intensity.  Generally  speaking,  a  higher  optical  intensity  generates  more  photo-carriers 
inside  the  PRC,  and  the  response  is  faster.31  Note  that  the  optical  power  density  used  here 
is  **20mW/mm2.  Also,  note  that  the  low  frequency  peaks  in  the  response  spectrum  are 
due  to  the  output  intensity  fluctuations  caused  by  the  source,  which  however  do  not  cause 
the  demodulator  to  drift  from  quadrature,  and  are  not  of  consequence  here. 

In  principle,  there  is  no  upper  limit  arising  from  the  TWM  process  for  the 
detectable  frequency  range  of  the  FBG  spectral  shift.  The  upper  frequency  is  limited 
only  by  the  FBG  sensor  response30  and  the  electronics  bandwidth  of  the  photodetector 
circuit.  As  shown  in  Figure  9,  a  high  frequency  dynamic  strain  signal  (a  580  kHz  3ps 
strain  applied  to  the  FBG  via  a  PZT-stretcher)  is  successfully  demodulated  by  the  TWM 
wavelength  demodulator. 
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Figure  9  A  580  kHz  3ps  strain  signal  demodulated  by  the  TWM  wavelength 

demodulator. 

Next  we  proceed  to  demonstrate  that  the  system  can  also  adapt  to  large 
temperature  drifts.  Temperature  drift  is  introduced  to  the  FBG  sensor  by  using  a  thermo¬ 
electric  cooler  (TEC)  module  connected  to  a  temperature  controller,  and  the  sensor 
temperature  is  monitored  by  a  thermistor.  A  0.1  Hz  sine  wave  is  used  as  input  to  the 
temperature  controller  and  the  temperature  cycle  amplitude  is  controlled  by  the  sine  wave 
amplitude.  The  temperature  drift  introduced  here  is  a  10°C  drift  within  5  seconds  (drift 
rate  of  2°C/s  which  is  extremely  high).  In  order  to  verify  that  the  temperature  drift  is 
effectively  introduced  to  the  FBG  sensor,  an  ANDO  (model  6317)  optical  spectrum 
analyzer  is  used  to  record  the  reflection  spectrum  of  the  FBG  sensor  as  the  temperature 
drift  is  introduced.  As  indicated  in  Figure  10,  the  temperature  drift  of  10°C  causes  a 
wavelength  shift  of  1 10pm,  which  is  close  to  the  expected  shift  of  130pm  (13pm/°C). 


Wavelength(nm) 

Figure  10  Reflection  spectra  of  the  FBG  sensor  subject  to  10°C  temperature  drift.  The 
center  wavelength  is  shifted  by  1 10pm  due  to  the  temperature  drift. 
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Figure  1 1  Time  trace  of  the  demodulator  amplitude  shows  no  variation  due  to  the  0.1  Hz 
10°C  temperature  cycling.  The  dashed  curve  is  what  the  demodulator  amplitude  would 
look  like  if  the  PRC  did  not  adapt  to  the  2.3  radian  phase  shift  caused  by  the  temperature 
drift.  Also  shown  is  the  thermistor  resistance  which  is  a  direct  measure  of  the  temperature 
change  in  the  FBG. 


Figure  1 1  is  the  time  trace  of  the  TWM  demodulated  signal,  which  shows  no 
variation  due  to  the  0.1  Hz  10°C  temperature  drift.  When  there  is  no  DC  field  applied 
onto  the  InP:Fe  PRC,  the  TWM  demodulator  signal  amplitude  was  originally  1.2  volts. 
This  increases  to  1.5  volts  due  to  the  diffraction  of  pump  beam  when  the  DC  field  is 
applied.  A  1 10pm  wavelength  shift  will  cause  an  equivalent  phase  shift  of  2.3  radians  for 
this  configuration.  If  such  a  large  phase  shift  were  to  be  picked  up  by  the  system,  the 
signal  beam  amplitude  would  undergo  a  large  amplitude  variation  as  indicated  by  the 
dashed  curve  in  Figure  11.  The  fact  that  no  such  shift  is  observed  indicates  that  the  TWM 
demodulator  adaptively  compensates  for  the  temperature  drift.  In  practical  applications, 
typical  environmental  temperature  drifts  will  be  on  the  order  of  a  few  degrees  per  hour, 
which  is  much  slower  than  the  case  considered  here  (2°C/s).  This  demonstrates  that  the 
TWM  wavelength  demodulation  system  will  be  able  to  adapt  to  all  reasonable 
environmental  temperature  drifts  of  FBG  sensors. 

2.6  Detection  of  transient  impact  signals  Next  we  demonstrate  the  ability  of  the  TWM 

wavelength  demodulator  to  monitor  transient 
(non-recurring)  events  such  as  acoustic  emissions,  vibrations  and  impact  signals.  This  is 
possible  becahse  the  TWM  demodulator  is  always  active  (unlike  scanned  demodulation 
schemes).  In  this  experiment,  an  FBG  sensor  was  covered  by  a  1mm  thick  diy  couplant 
made  from  silicon  epoxy.  A  3mm  thick  metal  plate  was  placed  on  top  of  the  couplant.  A 
3mm  diameter  ball  bearing  was  dropped  from  about  5  cm  above  the  metal  plate.  The 
oscilloscope  used  to  capture  the  demodulated  signals  was  triggered  by  the  demodulated 
signal  itself  with  a  trigger  level  set  to  be  slightly  higher  than  the  noise  level.  Figure  12 
shows  the  time  trace  of  the  impact  signal.  The  demodulated  FBG  sensor  signal  clearly 
shows  multiple  bouncing  of  the  ball  bearing.  The  frequency  of  the  first  impact  signal  is 
about  5  kHz. 
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Figure  12  Time  trace  of  the  impact  signal  induced  by  dropping  a  ball  bearing  on  a  metal 
plate.  The  inset  shows  the  experimental  configuration  of  the  impact  event. 


2.7  Four-channel  Multiplexed  TWM  Demodulator  addition  to  adaptivity  to  low 

frequency  drift,  another  major 
advantage  of  using  TWM  as  wavelength  demodulator  for  FBG  sensors  is  that  it  can  be 
multiplexed  without  significant  increase  in  cost.  This  is  because  all  the  channels 
(wavelengths)  share  the  same  PRC  and  there  is  no  expensive  feedback  electronics 
involved.  Multiplexed  TWM  systems  have  been  demonstrated  previously  at  NU  for 
phase-demodulation  applications.26"28  For  FBG  sensor  demodulation,  we  adopt  a  simpler 
wavelength  division  multiplexing  approach,  in  which  all  the  channels  share  the  same 
optical  setup  up  to  and  including  the  PRC,  and  subsequent  to  wave  mixing  in  the  PRC 
these  channels  are  separated  by  use  of  band-drop  filters  widely  used  in  the  optical 
communication  industry. 

Consider  multiple  FBG  sensors  with  distinct  spectral-reflectivities,  with  center 
wavelength  separation  AXC  that  is  chosen  to  be  sufficiently  large  so  that  stationary  optical 
interference  between  the  multiple  channels  cannot  occur.  In  this  case,  inside  the  PRC 
each  channel  creates  its  own  index  grating  of  different  grating  pitches.  The  change  in  the 
index  grating  pitches  is  related  to  the  channel  separation  and  the  signal  and  pump  beam 
angles,  and  this  can  be  shown  to  be: 


AA  = 


2sin(0/2)’ 


(9) 


in  which  AA  is  the  change  in  index  grating  pitch,  Alc  is  the  channel  separation  and  6  is 
the  angle  between  the  signal  and  pump  beams.  The  FBG  channel  spacing  Alc  should  be 
neither  too  small  to  avoid  crosstalk  involved  due  to  closely  packed  index  gratings,  nor  too 


large  to  maximize  the  number  of  sensors  that  can  be  used  within  the  limited  bandwidth  of 
the  C-band  (1530-1570nm).  In  this  work,  the  channel  separation  Ale  is  chosen  to  be 
4nm,  which,  according  to  Eq.  (9),  will  give  rise  to  an  index  grating  pitch  shift  of  76nm  at 
a  beam  angle  of  3°.  In  the  context  of  holographic  storage  elements  using  PRC’s,  it  has 
been  shown  that  multiple  gratings  can  be  written  in  a  single  PRC  with  negligible  cross¬ 
talk  if  the  grating  pitches  were  to  differ  by  as  little  as  0.03nm.32  We  therefore  expect  veiy 
small  crosstalk  in  our  setup,  if  any.  Note  that  a  4nm  channel  spacing  allows  10  channels 
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) 


in  the  C-band.  In  practice,  this  can  be  increased  by  using  narrower  channel  spacing,  at 
the  expense  of  decreased  dynamic  range. 


1548nm  1552nm  1556nm  1560nm 


Figure  13  Experimental  configuration  of  the  four-channel  two-wave  mixing  (TWM) 
wavelength  demodulator. 
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Figure  14  Simultaneous  demodulation  of  the  signals  from  four-FBG  sensors  using  a  4- 
channel  TWM  wavelength  demodulator.  FBG  sensor  1  (1548mn)  is  excited  with  1  kHz  5 
ps  strain,  sensor  2  (1552nm)  5  kHz  10  ps  strain,  sensor  3  (1556nm)  10  kHz  5  ps  strain 
and  sensor  4  (1550nm)  20  kHz  5  ps  strain. 
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As  a  proof  of  concept,  we  now  demonstrate  a  four-channel  TWM  wavelength 
demodulator.  The  experimental  configuration  is  shown  in  Figure  13.  Four  O.lnm  line- 
width  FBG  sensors  are  connected  in  series  and  are  centered  at  1548nm,  1552nm,  1556nm 
and  1560nm  respectively.  The  experimental  configuration  is  similar  to  that  of  the  single 
channel  configuration  shown  in  Figure  3  except  that  after  the  PRC,  there  is  a  free  space  to 
fiber  coupler  to  couple  the  free  space  light  into  a  set  of  four  band  drop  filters.  A  band 
drop  filter  is  a  three-port  device  that  transmits  a  certain  band  (i.e.  from  1546.8  nm  to 
1549.4nm)  and  reflects  all  the  other  wavelengths.  The  band  drop  filters  therefore 
decouple  the  TWM  demodulator  signals  from  the  various  FBG  sensors  prior  to 
photodetection.  The  width  of  the  banddrop  filter  should  be  chosen  to  be  wider  than  the 
expected  quasistatic  phase-shifts  of  the  FBG  sensors,  and  to  be  slightly  narrower  than  the 
FBG  channel  spacing. 

In  order  to  establish  the  crosstalk  of  the  four-channel  system,  we  applied  1 0  kHz  5 
ps  strain  on  FBG  sensor  1  (1548nm),  5  kHz  5  pe  on  FBG  sensor  2  (1552nm),  2  kHz  5  pe 
on  FBp  sensor  3  (1556nm)  and  20  kHz  5  pe  on  FBG  sensor  4  (1560nm)  simultaneously. 
Figure  14  shows  that  the  four  channels  can  be  demodulated  simultaneously.  Note  that  the 
low  frequency  fluctuations  seen  in  the  signals  are  due  to  environmental  noise-induced 
intensity  fluctuations  as  explained  earlier  in  this  paper,  and  these  can  be  removed  using  a 
balanced  detection  scheme  if  necessary.  Also  note  that  although  the  applied  signal 
amplitude  for  all  four  channels  are  the  same,  the  demodulated  signal  amplitude  for  each 
channel  is  different.  This  is  mainly  due  to  the  fact  that  each  channel  has  different  optical 
intensity  (due  to  non-uniform  EDFA  gain).  In  practice,  this  can  either  be  pre-calibrated, 
or  a  gain-flattened  EDFA  can  be  used. 


Frequency  (Hz) 

Figure  15  Fourier  spectrum  of  the  four-channel  TWM  wavelength  demodulated  signal 
shown  in  Figure  14.  There  is  no  detectable  crosstalk  between  these  four  channels. 

The  crosstalk  between  these  four  channels  can  be  inferred  from  the  Fourier 
spectra  of  the  signals  shown  in  Figure  14.  If  there  were  crosstalk,  we  would  expect  all  the 
frequency  components,  namely  2  kHz,  5  kHz,  10  kHz  and  20  kHz,  to  show  up  in  the 
signal  spectrum  of  each  channel.  Figure  15  is  the  Fourier  spectrum  of  the  four  channel 
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signals.  In  the  spectrum  of  channel  1,  there  is  only  the  expected  10  kHz  component  and 
no  other  frequency  components  are  observed.  This  is  true  also  for  the  other  three 
channels.  It  is  safe  to  conclude  that  the  crosstalk  between  these  four  channels  is  at  most 
comparable  to  the  noise  level  in  the  Fourier  spectrum,  which  is  at  least  30dB  below  the 
signal  level  in  our  case. 
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